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I.  INTRODUCTION 


Various  solid-fuel,  tubular  ramjet  projectiles  have  been  developed  at  the  U.S.  Army 
Ballistic  Research  Laboratory  (BRL).1,2  One  of  these  projectiles  has  been  designed  for 
use  as  a  tank  gun  training  round  (TGTR)  for  the  105mm,  M68  tank  cannon.  The  goal 
of  the  program  is  to  demonstrate  a  TGTR  of  low  dispersion  at  three  kilometers  and  a 
maximum  (safety)  range  of  eight  kilometers.  The  concept  of  the  TGTR  is  to  use  the 
thrust  (hence  low  drag)  of  the  solid-fuel  ramjet  (SFRJ)  projectile  to  obtain  a  ballistic 
match  with  low  drag  kinetic  energy  projectiles  up  to  three  kilometers.  Upon  depletion 
of  the  solid  propellant  and  choking  of  the  internal  flow,  the  SFRJ  will  become  a  high- 
drag  projectile  with  limited  range.  The  successful  demonstration  of  the  SFRJ  projectile 
has  included  the  verification  of  flight  and  motor  performance,  limited  safety  range  and 
ballistic  match  to  a  kinetic  energy  projectile  up  to  2.5  kilometers  with  somewhat  higher 
dispersion  than  desired.  Investigations  into  the  combustion  efficiency  are  currently  under 
way  in  an  effort  to  extend  the  ballistic  match  to  3  kilometers.  Computational  flow  modeling 
significantly  contributes  to  this  investigation  by  adding  to  the  current  understanding  of 
the  SFRJ  internal  flow. 

Computational  modeling  of  the  internal  and  external  flow  for  a  75mm  SFRJ  projectile 
has  been  under  way  at  the  BRL  Launch  and  Flight  Division  in  recent  years.  A  computa¬ 
tional  fluid  dynamics  (CFD)  code  utilizing  an  implicit,  factored,  time-stepping  algorithm 
in  a  zonal  grid  framework  has  been  developed  by  Chakravarthy.3  This  code  employs  a 
class  of  numerical  algorithms,  termed  total  variational  diminishing  or  TVD,  which  do  not 
require  the  inclusion  of  smoothing  or  dissipation  functions  to  achieve  numerical  stability. 
The  code  can  be  used  in  conjunction  with  various  turbulence  and  separated  flow  modeling 
techniques.  Modeling  strategies  investigated  include:  a)  the  Baldwin-Lomax  turbulence 
model  4  applied  throughout  the  flowfield;  b)  the  Baldwin-Lomax  model  applied  outside  of 
backflow  regions,  and  the  backflow  turbulence  model  of  Goldberg  5  applied  within  these 
regions.  This  code  has  been  previously  employed  in  the  solution  of  subsonic,  transonic, 
supersonic  and  mixed  flow  problems  including  complex  supersonic  inlet  and  nozzle  flows 

1  Mermagen  W.H.,  )  a  la  man  chi  L  R.J. ,  "First  Diagnostic  Tests  of  a  75mm  Solid  Fuel  Ramjet  Tubular  Projectile,''  US  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-03263,  June  19SS.  (AD  A1S0596) 

2 Mermagen  W.H.,  Yalamanckih  R.J.,  "Experimental  Tests  of  a  t05/75mm  Solid  Fuel  Ramjet  Tubular  Projectile,"  US  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-MR-S416,  December  1964.  (AD  B069766) 

3  Chakravarthy,  S.R.  "A  New  Computational  Capability  for  Ramjet  Projectiles,"  ARBRL-CR-595,  U.S.  Army  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  March  1966  (Also  Goldberg,  U.,  Chakravarthy,  S.,  and  Nus ce,  M., 
AIAA-67-241 1,  Proceedings  of  the  14th  AIAA  Atmospheric  Flight  Mechanics  Conference,  Monterrey,  CA,  August  17-19, 
1967.) 

4  Baldwin,  B.S.  and  Lomax,  H.,  “ Thin  Layer  Approximation  and  Algebraic  Model  for  Separated  Turbulent  Flows,”  AIAA- 
76-257,  Proceedings  of  the  15th  AIAA  Aerospace  Sciences  Meeting,  Huntsville,  AL,  January  15-16,  1976. 

s  Goldberg  U.C.,  “Separated  Flow  Treatment  with  a  New  Turbulence  Model,"  AIAA  Journal.  Vo/.  24,  No.  10,  October 
1966,  pp.  1711-1713. 
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by  Chakravarthy  et.  al.  6-11  This  report  describes  application  of  the  code  to  an  inert 
SFRJ  projectile.  Results  for  axisymmetric  (zero  yaw)  cold  flow  (no  combustion  or  mass 
injection)  are  presented. 

Previous  computational  efforts  include  the  application  of  an  unsteady,  thin-layer 
Navier-Stokes  code  to  ramjet  type  configurations12  and  an  axisymmetric  method  of  char¬ 
acteristics  solution  for  ramjet  inlet  flow.  13  The  latter  solution  assumed  separated  flow  in 
the  inlet  and  an  experimentally  determined  inlet  pressure  level. 

The  SFRJ  flight  vehicle  had  a  1.7  inch  diameter  injector  and  a  1.54  inch  diameter 
nozzle.  During  wind  tunnel  tests,  small  changes  in  total  pressure  produced  dramatically  dif¬ 
ferent  internal  wall  pressure  distributions  that  were  not  observed  for  other  injector/nozzle 
diameter  combinations.  14’15  This  unusual  behavior  could  reflect  bistable  flow  or  highly 
complex  flow  states  for  this  configuration.  Comparisons  between  computed  and  measured 
internal  wall  pressures  jure  made  for  several  injector/nozzle  diameter  combinations  and  a 
single  total  pressure  for  which  stable  wall  static  pressures  were  observed.  The  detailed  anal¬ 
ysis  of  the  effect  of  different  total  pressures  for  the  1.7/1.54  combination  is  not  addressed 
in  this  report. 

Cold  flow,  zero-yaw  computations,  for  which  internal  flow  measurements  exist,  are  a 
necessary  prelude  to  further  extensions  of  the  computational  analysis  for  the  SFRJ  pro¬ 
jectile.  The  zonal  griding  approach  of  the  code  allows  incorporation  of  a  combustion 
model,  such  as  that  described  by  Vos,16  into  the  combustion  section  of  the  internal  geom¬ 
etry.  After  testing  and  validation,  these  reacting  flow  computations  can  be  extended  to  the 
three-dimensional  case,  non-zero  yaw,  with  direct  comparisons  with  measured  aerodynamic 
drag. 

®  Chakravarthy  S.R.,  Szema  A.}'.,  Goliber g  U.C.,  Gonki  J.J.  (Rockwell  International  Science  Center)  and  Other  S.  (Uni- 
versity  of  California),  “Application  of  a  New  Class  of  High  Accuracy  TVD  Schemes  to  the  Navier-Stokes  Equations AIAA- 

85- 0165,  Proceedings  of  the  23rd  AIAA  Aerospace  Sciences  Meeting,  Reno  NV.,  January  14-17,  1985. 

7  Chakravarthy  S.R.,  Szema  K.Y.,  “An  Euler  Solver  for  Three-Dimensional  Supersonic  Flows  with  Subsonic  Pockets ,” 
AlAA-65-1 70S,  Proceedings  of  the  18th  AIAA  Fluid  Dynamics  and  Plasmadynamics  and  Lasers  Conference ,  Cincinnati  OH ., 
July  16-18,  1985. 

*  Chakravarthy  S.R.,  “The  Versatility  and  Reliability  of  Euler  Solvers  Based  on  High- Accuracy  TVD  Formulations,  ”  A1AA- 

86- 0243,  Proceedings  of  the  24th  AIAA  Aerospace  Sciences  Meeting,  Reno  NV.,  January  6-9,  1986. 

9  Chakravarthy  S.R.,  Harten  A.,  Osher  S.,  “ Essentially  A  on- Oscillatory  Shock-  Capturing  Schemes  of  Arbitrarily- High  Ac¬ 
curacy,”  AIA A-86-0339,  Proceedings  of  the  24th  AIAA  Aerospace  Sciences  Meeting,  Reno  NV.,  January  6-9,  1986. 

10 Sova  G.J.,  “Application  of  I^aner-Stokes  Analysis  to  Predict  the  Internal  Performance  of  Thrust  Vectoring  Two- 
Dimensional  Convergent- Divergent  A ozzlts,”  A 1AA-88-2S86,  Proceedings  of  the  6ih  AIAA  Applied  Aerodynamics  Conference, 
Williamsburg,  \'A ,  June  6-8,  1588. 

11  Chakravarthy,  S.R.,  I\.-Y.  Szema,  and  Haney,  J.W.,  “Unified  “Nose-to-Tail”  Computational  Method  for  Hypersonic 
Vehicle  Applications ,”  AlAA-66-2564,  Proceedings  of  the  6th  AIAA  Applied  Aerodynamics  Conference,  Williamsburg,  V A , 
June  6-8,  1988. 

** Kietubicz,  C.J..  and  Heavey,  h.R.,  “Computational  Flow  Field  Predictions  for  Ramjet  and  Tubular  Projectiles,”  Proceed¬ 
ings  of  the  8th  ADPA  International  Symposium  on  Ballistics,  Orlando,  FL,  October  23-25,  1964. 

^Danfcerp,  J.E.,  and  Sigal,  A.,  “Evaluation  of  Solid  Fuel  Ramjet  Projectile  Aerodynamic  Characteristics,”  Proceedings  of 
the  10th  ADPA  International  Symposium  on  Ballistics,  San  Diego,  CA,  October  27-29,  1967. 

14  Kayser  L.D.,  Yalamanchilt  R.J.,  Trexler  C.,  “Pressure  Measurements  on  the  Interior  Surface  of  a  75mm  Tubular  Projectile 
at  Mach  4t  VS  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland,  report  in  preparation. 

15  Yalamanchili  R.J.,  unpublished  wind  tunnel  data  for  the  75mm  Tubular  Projectile  at  Mach  4,  US  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground ,  Maryland. 

16  Vos,  J.B.,  “Calculating  Turbulent  Reacting  Flows  Using  Finite  Chemical  Kinetics,”  AIAA  Journal.  Vol.  25,  No.  10, 
October  1967,  pp. 1365- 1372. 
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II.  COMPUTATIONAL  APPROACH 


The  present  CFD  approach  can  be  used  to  predict  the  compressible  flowfield  in  and 
around  an  aerodynamic  projectile  by  solving  the  3D  Navier-Stokes  equations.  At  present 
only  axisymmetric  (zero  yaw)  flows  have  been  investigated.  These  equations  are  solved 
with  the  assumption  that  the  flow  medium  is  air  behaving  as  a  perfect  gas  and  that  no 
chemical  reactions  are  occurring.  As  a  result  only  inert  SFRJ  projectiles  are  addressed  at 
present.  Both  laminar  and  turbulent  flows  are  investigated,  thus  an  adequate  turbulence 
model  is  required  for  closure.  In  addition,  backflow  regions  can  be  present,  thus  a  back- 
flow  turbulence  model  is  included.  The  equations  are  transformed  into  conservation  law 
form  and  discretized  using  finite  volume  approximations  and  the  TVD  formulation.  The 
resulting  set  of  equations  is  solved  using  an  implicit,  factored,  time-stepping  algorithm. 
This  solution  takes  place  on  a  computational  grid  that  is  generated  around  the  projectile 
in  zones  where  the  zonal  boundaries  can  be  made  transparent  to  the  flowfield  calculation. 


1.  EQUATIONS  OF  MOTION. 


The  compressible,  Reynolds- averaged  Navier-Stokes  equations  for  2D / axisymmetric 
flow  are  written  in  the  following  conservation  form.  The  dependent  variables  u,  v  and  e 
are  mass- averaged,  with  e  being  the  specific  total  internal  energy,  T  being  temperature,  p 
and  p  being  the  mean  density  and  pressure,  respectively,  and  t  being  time. 
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where  a  =  1  for  axisymmetric  flow  and  0  for  two  dimensional  flow. 


(7) 

(8) 
(9) 


In  Equations  2-7,  the  laminar  and  eddy  viscosities,  p  and  pt,  are  implicitly  divided  by 
the  reference  Reynolds  number.  The  equations  used  for  the  Euler  (inviscid)  calculations 
are  obtained  from  Equations  1-9  by  setting  both  laminar  and  eddy  viscosities  to  zero.  In 
all  calculations,  the  flow  medium  (air)  was  assumed  to  be  a  perfect  gas,  satisfying  the 
equation  of  state 

p  =  pRT  (10) 

The  following  power  law  was  used  to  relate  molecular  viscosity  to  temperature:  17 


h_ 

ho 


(11) 


where  pa  =  0.1716  mP,  T0  —  491.6  R,  and  n  =  0.64874.  The  laminar  and  turbulent  Prandtl 
numbers,  Pr  and  Prt,  were  assumed  constant  with  values  of  0.72  and  0.9  respectively.  The 
ratio  of  specific  heats,  7,  was  also  assumed  constant  and  equal  to  1.4.  The  specific  heat 
capacities  at  constant  volume  and  pressure,  Cv  and  Cp,  are  related  as  7  =  Cv/Cv. 


In  the  £  —  r)  computational  plane,  Equation  1  is  transformed  into  the  finite  volume 
conservation  law  form  represented  by  6 

fr  +  ((*f  -  X-G,< +(-*F + X<G)” +  7  -  7)  = 0  (12) 

where  £  and  7  are  the  new  independent  variables  and  xv,  y^,  and  yv  are  the  four 
transformation  coefficients  obtained  numerically  from  the  mapping  procedure.  The  “Area” 
in  Equation  12  denotes  the  area  of  the  finite  volume  cell  under  consideration  at  the  time 
of  discretization  of  the  equations.  The  transformed  time  variable  is  represented  by  r. 


2.  TURBULENCE  MODELING. 

The  internal  flowfield  of  the  SFRJ  projectile  cam  include  large  regions  of  recirculatory 
flow,  induced  by  both  shock  waves  and  by  sharp  geometrical  discontinuities.  Indeed  in¬ 
ternal  surface  pressure  measurements  made  on  an  instrumented  SFRJ  wind  tunnel  model 
(unfueled)  indicate  that  such  regions  do  exist  in  the  inlet  and  combustion  sections.  1415 
Good  modeling  of  these  regions  is  critical  to  the  overall  flowfield  solution  quality.  However 
most  existing  turbulence  models  either  do  not  treat  such  regions  or  do  so  in  an  ad  hoc 
fashion  that  is  frequently  inadequate.  A  notably  different  approach  is  the  use  of  a  full 
Reynolds  stress  closure  model,  involving  the  solution  of  five  coupled  partial  differential 

17 Motor  G.,  Ben-Dor  G.,  and  Igra  0.,  “A  Simple  and  Accurate  Erprettion  for  the  Viacoiity  of  Nonpolar  Diatomic  Gael 
%p  to  10,000  K.  ”  AIAA  Journal.  Vol.  23,  A :o.  4,  April  1985,  pp.  636-638. 
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equations  (for  two-dimensional  flows)  for  the  three  normal  stresses,  the  shear  stress,  and 
the  length  scale.  These  must  be  supplemented  by  a  wall  function  to  provide  turbulence 
quantities  across  viscous  regions  adjacent  to  solid  surfaces.  Such  a  wall  function  is  usu¬ 
ally  some  form  of  the  law-of-the-wall,  which,  according  to  experimental  observations,  does 
not  apply  to  detached  flows.  Thus,  an  expensive  and  time  consuming  computation  of  the 
Reynolds  stresses  is  coupled  with  a  questionable  near-wall  formulation. 

To  improve  the  predictive  capability  of  separated  flows  using  current  Reynolds-averaged 
Navier-Stokes  codes,  a  new  turbulence  model  has  been  recently  developed.5  The  new  tur¬ 
bulence  model  is  based  on  experimental  observations  of  detached  flows.  The  model  pre¬ 
scribes  turbulence  kinetic  energy  ( k )  and  its  dissipation  (e)  analytically  within  backflows. 
A  Gaussian  variation  of  k  normal  to  walls  is  assumed.  The  length  scale  of  turbulence 
is  proportional  to  the  local  distance  from  the  wall  to  the  edge  of  the  viscous  sublayer, 
which  is  located  outside  the  backflow  region.  The  latter  feature  is  a  basic  assumption  of 
the  model.  The  stress  scale  is  the  local  maximum  Reynolds  stress,  which  typically  occurs 
around  the  middle  of  the  boundary  layer,  well  outside  the  separation  bubble.  This  scale 
must  be  supplied  by  a  turbulence  model  that  is  used  beyond  backflow  regions. 

The  main  equations  of  the  backflow  model  are  given  in  Reference  5.  A  formula  for  the 
eddy  viscosity  distribution  within  backflows  results  and  is  used  to  supply  eddy  viscosity  for 
the  Reynolds-averaged  equations  when  the  calculations  are  done  inside  separation  bubbles. 
Outside  of  them,  another  turbulence  model  (for  example  Baldwin- Lomax4)  supplies  the 
values  of  eddy  viscosity.  While  the  Baldwin-Lomax  turbulence  model  is  used  to  detect 
flow  separation  and  to  initiate  application  of  the  backflow  model,  the  latter  model  may 
relocate  the  separation  point  (and  the  reattachment  point,  if  one  exists).  The  location 
of  flow  separation  as  predicted  by  the  Baldwin-Lomax  turbulence  model  and  Goldberg’s 
backflow  model  are  significantly  different  for  the  applications  reported  in  Reference  18. 
For  further  details  of  how  the  model  treats  the  influence  of  large  eddies  residing  outside 
detached  regions,  the  history  effect  of  these  eddies  downstream  of  reattachment,  and  the 
mutual  influence  of  multiple  walls  on  the  eddy  viscosity,  see  Reference  18.  The  model 
has  been  tested  successfully  for  a  variety  of  flow  conditions  and  body  geometries  including 
a  backward-facing  step  in  subsonic  flow,18  an  axisymmetric  bump  in  transonic  flow,18  and 
an  axisymmetric  boattail  in  transonic  flow.19 


3.  COMPUTATIONAL  ALGORITHM. 


The  spacial  discretization  technique  for  the  equations  of  motion  must  be  reliable  and 
robust  if  it  is  to  successfully  capture  the  complex  physics  of  SFRJ  internal  and  external 
flowfields.  The  TVD  formulation  for  the  convection  terms  (the  hyperbolic  part  of  the 
time-dependent  Navier-Stokes  equations),  along  with  a  special  treatment  of  the  diffusion 
terms,  provides  an  appropriate  simulation.  In  order  to  enable  the  solution  procedure,  a  time 
discretization  operator  has  to  be  utilized  as  well.  Any  conventional  method,  suitable  for  the 

14 Goldberg  U.C.,  “Prediction  of  Separated  Flow  WilA  A  New  Turbulence  Model,”  Proceeding/  of  the  Sth  International 
Conference  on  Numerical  Method »  in  Laminar  and  Turbulent  Flow.  Montreal,  Quebec,  Canada,  July  6-10,  1987.  (See  also 
A1AA  Journal.  Vol.  t6.  No.  d,  April  1988,  pp.  405-408). 

18  Goldberg  U.C.,  “Separated  Flow  Calculation/  With  A  New  Turbulence  Model,”  Proceeding/  of  the  Fir/t  World  Congre// 
on  Computational  Mechanics  the  University  of  Texas  at  Austin,  September  22-26,  1986. 
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Navier-StOKes  equations,  can  be  used  together  with  the  space  discretization  methodology 
described  above.  This  includes  approximate  factorization  and  relaxation  techniques.  In 
recent  years,  TVD  formulations  have  been  constructed  for  shock-capturing  finite- difference 
methods.6"9  Near  large  gradients  in  the  solution  (extrema),  TVD  schemes  automatically 
reduce  to  first-order  accurate  discretizations  locally,  while  away  from  extrema  they  can  be 
constructed  to  be  of  higher-order  accuracy.  This  local  effect,  which  is  necessary  to  prevent 
the  total  variation  from  increasing,  restricts  the  maximum  global  accuracy  possible  for 
TVD  schemes  to  third  order  for  steady-state  solutions. 

These  methods  manifest  many  properties  desirable  in  numerical  solution  procedures. 
By  design,  they  avoid  numerical  oscillations  and  “expansion  shocks”  while  at  the  same  time 
being  higher-order  (more  than  first-prder)  accurate.  (“Expansion  shocks”  are  shock  waves 
which  do  not  satisfy  the  entropy  inequality).  TVD  formulations  are  also  based  on  the 
principle  of  discrete  or  numerical  conservation,  which  is  the  numerical  analog  of  physical 
conservation  of  mass,  momentum,  and  energy.  This  results  in  TVD  schemes  being  able  to 
“capture”  discontinuities  with  ease  and  high  resolution.  At  a  fundamental  level,  they  are 
based  on  upwind  schemes;  therefore,  they  closely  simulate  the  signal  propagation  properties 
of  hyperbolic  equations. 

Schemes  based  on  the  TVD  formulation  are  completely  defined.  In  other  words,  the 
user  does  not  have  to  specify  any  “numerical”  parameters  such  as  “dissipation  coefficients” 
to  be  able  to  carry  out  the  computations.  In  contrast,  central  difference  schemes  involve 
dissipation  terms  for  stability  and  have  one  or  more  coefficients  that  must  be  judiciously 
chosen  to  achieve  desirable  results.  Methods  based  on  TVD  formulations  represent  mature 
CFD  technology  and  are  being  applied  to  a  wide  class  of  problems.  8 

Proper  treatment  of  the  diffusion  terms  of  the  Navier-Stokes  equations  is  also  impor¬ 
tant  in  the  construction  of  reliable  numerical  methods.  Unidirectional  second  derivative 
terms  are  treated  by  using  central  difference  approximations.  Cross  derivatives  are  repre¬ 
sented  by  finite-differences,  the  nature  of  which  depends  upon  the  sign  of  the  coefficient  of 
such  terms.  This  treatment  augments  diagonal  dominance  of  the  resulting  set  of  discretized 
equations,  without  detracting  from  the  accuracy  and  while  adding  to  the  reliability  of  the 
numerical  procedure.  Further  details  can  be  found  in  Reference  6.  The  diagonal  domi¬ 
nance  of  TVD  schemes  (with  suitable  discretized  cross  derivative  terms)  makes  it  possible 
to  use  relaxation  methods.  However,  conventional  approximate  factorization  methods  are 
used  to  solve  multidimensional  implicit  formulations  in  the  present  study. 


4.  COMPUTATIONAL  GRID. 

The  problem  of  computing  internal/external  SFRJ  projectile  flowfields  is  complicated 
by  the  internal  geometry  involved.  As  illustrated  in  Figure  1,  the  SFRJ  consists  of  sev¬ 
eral  sharp  corners  that  would  severely  hamper  conventional  grid  generation  schemes  that 
require  one  set  of  grid  lines  to  be  tangent  to  the  surface  and  another  set  to  be  normal 
to  it.12  The  SFRJ  geometry  is  more  easily  gridded  by  the  zonal  approach.  The  internal 
geometry  of  the  SFRJ  in  broken  up  into  four  zones  of  simple  geometric  shape  (Figure  2). 
In  each  zone  any  approach  to  generating  the  grid  can  be  used,  such  as  algebraic  methods, 
differential  equations  methods,  etc.  In  this  case  a  simple  algebraic  grid  is  used  with  grid 
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clustering  near  the  surface  and  other  regions  where  high  gradients  are  expected  (Figure 
3).  In  the  zonal  approach,  the  computational  method  and  computer  program  are  con¬ 
structed  in  such  a  manner  that  each  zone  may  be  considered  as  am  independent  module, 
interacting  with  each  other  before  or  after  the  information  corresponding  to  each  zone  is 
updated  one  cycle.  In  addition,  the  zonal  boundaries  cam  be  made  transparent  to  flowfield 
phenomena  (e.g.  shock  waves).  The  actual  grid  used  for  these  computations  consisted  of 
the  following  dimensions  for  zones  1  through  6:  60x61,  7x30,  80x49,  40x30,  15x49,  amd 
45x20.  These  grid  dimensions  represent  a  refinement  over  a  coarse  grid  that  was  initially 
used  (30x35,  7x17,  50x30,  30x17,  15x25,  45x20). 3  The  computationad  results  were  found  to 
be  essentiailly  independent  of  these  two  grids  when  all  other  factors  (e.g.  flow  turbulence 
and  transition  modelling)  were  the  same.  However,  the  fine  grid  wais  preferred  for  resolving 
flow  details  within  boundary  layers.  The  grids  have  been  generated  in  such  a  manner  that 
the  horizontad  faunilies  of  mesh  lines  within  the  projectile  are  continuous  between  zones 
and  the  verticad  faimilies  are  continuous  in  the  zones  external  to  the  SFRJ.  While  this  type 
of  grid  is  natural  for  this  problem,  the  computationad  methodology  permits  a  wider  class 
of  patched  zonad  grids  where  neither  family  of  grid  lines  need  be  continuous  across  zones. 

III.  RESULTS 

The  computational  solutions  are  compared  with  internal  surface  pressure  measure¬ 
ments  of  an  inert  SFRJ.  The  model  was  instrumented  with  pressure  taps  amd  mounted 
in  the  Mach  4,  nine-inch  blowdown  tunnel  at  the  NASA  Lamgley  Research  Center.  The 
details  of  these  tests  are  described  in  Reference  14.  The  freestream  Mach  number  and 
Reynolds  number  were  4.03  amd  about  20  million  per  foot,  respectively.  The  model  was  at 
zero  yaw. 

In  the  computational  solutions,  the  freestream  Mach  number  was  increased  gradually 
as  the  numerical  iterations  progressed.  At  a  given  point  in  the  iteration  cycle,  the  Mach 
number  was  fixed  and  the  iterations  continued  until  convergence.  For  the  Euler  (inviscid) 
solutions,  the  Mach  number  was  linearly  increased  from  0.03  to  4.03  for  3000  iterations  and 
then  held  fixed  for  the  balance  of  5000  total  iterations.  For  the  Navier-Stokes  solutions,  the 
Mach  number  was  linearly  increased  from  0.03  to  4.03  for  1000  iterations  and  then  held  fixed 
for  the  balance  of  5000  total  iterations.  The  more  conservative  treatment  for  the  Euler  cases 
was  necessary  given  the  complex  geometry  of  the  SFRJ  and  the  lack  of  natural  viscosity  to 
damp  transients  developing  in  the  inviscid  flow.  In  all  cases  convergence  was  determined 
by  the  arrival  at  a  time- asymptotic  steady-state  solution.  The  CFL  number  (stability 
requirement  for  the  numerical  solution  of  hyperbolic  partial  differential  equations)  for 
these  results  was  set  at  1.0  for  the  Euler  solutions,  and  between  .5  to  2.5  for  the  Navier- 
Stokes  solutions,  although  solutions  can  be  achieved  in  about  1000  iterations  with  CFL 
numbers  as  high  as  10.  Typical  run  times  on  a  Cray  XMP  4/8  computer  for  Euler  and 
Navier-Stokes  solutions  were  0.8  and  1.2  CPU  hours,  respectively. 

The  Navier-Stokes  solutions  consisted  of  both  laminar  and  turbulent  cases.  In  the 
turbulent  cases  the  Baldwin-Lomax  model  was  used  both  with  and  without  the  Goldberg 
backflow  model  previously  described.  For  those  cases  in  which  the  latter  model  was  in¬ 
voked,  it  was  only  used  if  the  computed  flow  conditions  (based  on  the  Baldwin-Lomax 
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model)  indicated  flow  separation.  If  this  condition  did  not  occur  the  Baldwin-Lomax 
model  was  used  exclusively.  Due  to  the  zonal  nature  of  the  code,  a  flow  type  (laminar  or 
turbulent)  can  be  assigned  to  each  zone.  In  certain  solutions  all  zones  were  solved  with 
the  same  flow  type  assumption. 

1.  IMPORTANCE  OF  FLOW  TRANSITION  AND  TURBULENCE. 

As  indicated  previously,  the  present  CFD  code  can  be  run  in  severed  modes.  These 
modes  include  inviscid,  laminar,  and  turbulent  flow  and  involve  the  solution  of  the  Euler  or 
the  Navier-Stokes  equations.  For  any  given  geometry  and  flow  conditions,  the  mechanisms 
that  control  transition  from  laminar  to  turbulent  flow  and  the  presence  of  attached  or 
detached  flow  are  not  completely  understood.  To  achieve  reliable  results  then,  the  code 
must  be  used  in  several  modes.  Comparison  of  the  results  can  indicate  the  impact  of 
turbulence  and  flow  detachment  for  a  particular  configuration.  In  addition,  various  flow 
transition  locations  can  be  studied. 

Figure  4  shows  computed  internal  wall  pressure  distributions  for  the  configuration 
with  1.7  inch  injector  and  1.6  inch  nozzle  diameters.  The  code  was  run  assuming  inviscid, 
laminar,  turbulent,  and  turbulent /detached  flow.  For  the  turbulent  cases,  transition  from 
laminar  to  turbulent  flow  was  specified  near  the  leading  edge  of  the  inlet.  The  best  overall 
agreement  of  the  computations  with  pressure  data  in  the  combustion  section  (3  <  x  <  9) 
is  achieved  when  the  backflow  turbulence  model  is  used.  In  this  section,  the  flow  can  be 
characterized  as  turbulent  and  detached.  Studies  of  the  computed  flow  velocity  vectors 
indicate  that  the  combustion  region  is  dominated  by  a  separation  bubble.  Here,  the  flow 
separation  point  is  determined  by  the  sharp  corner  of  the  injector. 

In  the  inlet  (0  <  x  <  3)  the  computations  with  turbulence  predict  an  expansion 
at  the  leading  edge,  followed  by  a  gradual  compression  that  starts  at  about  one-third 
the  inlet  length  from  the  leading  edge  (i.e.  the  flow  separation  point).  The  Baldwin- 
Lomax  turbulence  model  predicts  the  separation  point  slightly  further  upstream  than  the 
backflow  turbulence  model  of  Goldberg.  However,  the  wind  tunnel  data  show  a  constant 
pressure  over  the  entire  inlet  indicating  flow  separation  at  or  near  the  leading  edge  of  the 
inlet.  Unexpectedly,  the  computations  assuming  laminar  flow  better  match  the  data  in 
the  inlet.  This  result  suggests  either;  1)  that  the  inlet  flow  is  turbulent  and  detached 
but  the  Baldwin-Lomax  turbulence  model  and  the  Goldberg  backflow  model  predict  flow 
separation  further  downstream  than  indicated  by  the  wind  tunnel  data,  or  2)  that  the  inlet 
flow  is  transitional  (i.e.  neither  fully  turbulent  nor  laminar)  and  the  point  of  transition 
from  laminar  to  turbulent  flow  cannot  be  arbitrarily  specified. 

Figure  5  shows  the  effects  of  flow  transition  location  on  turbulent/detached  flow  com¬ 
putations  for  the  configuration  with  1.7  in.  injector  and  1.6  in.  nozzle  diameters.  For 
the  two  computations  shown,  flow  transition  has  been  specified  at  either  the  beginning 
of  the  inlet  or  at  the  end  of  the  inlet  (i.e.  at  the  injector).  For  the  latter  case,  the  inlet 
flow  is  computed  as  laminar  and  the  best  overall  comparison  with  wind  tunnel  data,  in 
this  region,  is  achieved.  Studies  of  the  computed  flow  velocity  vectors  indicate  that  a 
large  separation  bubble  extends  over  the  entire  inlet  length.  The  laminar  flow  assumption 
more  accurately  predicts  the  separation  point  in  the  inlet.  This  result  is  unexpected  since 


laminar/separated  flows  are  unusual.  However,  as  stated  above  the  turbulence  models  I 

employed  may  be  deficient  in  predicting  the  flow  separation  point  in  the  inlet.  Specifying  I 

transition  at  the  injector  improves  agreement  with  data  behind  the  injector  but  degrades  I 

agreement  at  the  end  of  the  combustion  section.  For  the  entire  configuration,  the  best  ^ 

overall  agreement  between  computed  and  measured  wall  pressures  is  achieved  for  injector 
transition. 

Figure  6  shows  the  results  of  a  similar  flow  transition  study  for  a  configuration  with 
1.9  in.  injector  and  1.1  in.  nozzle  diameters.  In  contrast  to  the  1.6  in.  nozzle,  the  smaller  j 

nozzle  diameter  causes  a  normal  shock  at  the  leading  edge  of  the  SFRJ  with  subsonic 
flow  throughout  the  interior.  For  this  configuration,  the  best  overall  agreement  with  wind 
tunnel  data  occurs  when  flow  transition  is  specified  at  the  beginning  of  the  inlet.  In 
this  case,  the  normal  shock  has  determined  the  point  of  flow  separation.  In  summary, 
the  SFRJ  configurations  studied  to  date  indicate  that  flow  transition  may  occur  at  the 
beginning  of  the  inlet  for  subsonic  cases  (nozzle  diameter  1.1  in.)  and  at  the  end  of  the  ^ 

inlet  for  supersonic  cases  (nozzle  diameter  1.6  in.).  ] 


2.  DISCUSSION  OF  COMPUTED  RESULTS. 

Figure  7  shows  comparisons  of  internal  surface  pressure  measurements  and  turbu¬ 
lent/detached  flow  computations  for  three  SFRJ  configurations.  The  configuration  used 
in  SFRJ  flight  tests  (1.7  inch  injector  and  1.54  inch  nozzle  diameter)  is  included.  These 
configurations  differ  only  in  the  nozzle  diameter;  1.6  in.,  1.54  in.,  and  1.1  in.  The  injector 
diameter  was  1.7  in.  in  each  case.  The  criterion  established  in  the  previous  section  was 
used  to  specify  the  flow  transition  location.  The  lower  pressure  level  for  the  1.6  and  1.54 
in.  nozzle  configurations  is  indicative  of  supersonic  internal  flow,  with  the  exception  of 
the  wall  boundary  layers  and  detached  regions.  Although  obscured  by  the  scale  of  this 
plot,  the  overall  agreement  between  computation  and  measurements  for  the  smaller  nozzle 
diameters  is  notably  better  than  that  for  the  1.1  in.  nozzle  (see  Figure  5).  Figure  S  shows 
the  pressure  contours  for  the  1.6  in.  nozzle  SFRJ.  An  oblique  shock  is  attached  to  the 
leading  edge  of  the  SFRJ  and  intersects  the  centerline  of  the  model  as  a  normal  shock. 
Examination  of  the  computed  velocity  vectors  shows  that  the  inlet  flow  is  dominated  by 
a  separation  bubble  that  terminates  at  the  injector  where  an  oblique  shock  is  generated. 
The  flow  in  the  combustion  section  is  also  dominated  by  a  separation  bubble  that  extends 
from  the  injector  to  the  nozzle  and  well  into  the  flowfield.  A  zero  wall  pressure  gradient 
extends  throughout  this  region  (see  Figure  8).  Several  shock  reflections  occur  throughout 
the  combustion  and  nozzle  sections. 

Danberg,13  used  an  axisymmetric  method  of  characteristics  solution  to  investigate 
possible  inlet  flow  patterns.  This  analysis  assumed  separated  flow  over  the  entire  length 
of  the  inlet  and  a  constant  inlet  wall  pressure  at  the  experimentally  determined  pressure 
level.  This  calculation  showed  that  an  oblique  shock  emanates  from  the  leading  edge  of 
the  inlet  and  terminates  in  a  normal  shock  at  the  centerline.  Danberg  showed  that  the 
inlet  flow  pattern  bears  a  striking  resemblance  to  that  of  an  over-expanded  jet  (where  the 
jet  exit  coincides  with  the  inlet  entrance). 

The  higher  pressure  level  for  the  1.1  in.  nozzle  diameter  (Figure  7)  indicates  subsonic 
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internal  flow  caused  by  an  expelled  leading  edge  normal  shock.  The  computed  wall  pressure 
matches  the  distribution  but  not  the  level  of  the  measured  pressure.  Indeed  the  computed 
pressure  level  better  matches  the  measured  result  for  the  1.2  in.  nozzle  (see  Figure  10, 
Reference  3).  This  would  indicate  a  nozzle  wall  boundary  layer  thickness  which  is  predicted 
to  be  thinner  that  the  actual  one.  This  could  be  the  result  of  turbulence  model  deficiency 
since  these  computed  pressures  were  found  to  be  independent  of  grid  resolution  normal  to 
the  nozzle  wall.  Investigations  with  other  turbulence  models  should  clarify  this  point. 

Figure  9  shows  the  pressure  contours  for  the  1.1  in.  nozzle  SFRJ.  A  normal  shock  is 
clearly  seen  at  the  leading  edge.  Upon  detailed  examination,  the  presence  of  flow  spillage 
to  the  exterior  has  been  found.  The  subsonic  interior  flow  expands  through  the  diverging 
inlet,  and  the  pressure  remains  relatively  constant  throughout  the  combustion  section. 
Expansion  back  to  supersonic  flow  is  accomplished  by  the  nozzle.  Examination  of  the 
computed  velocity  vectors  shows  that  the  inlet  flow  separates  at  about  two-thirds  of  the 
length  from  the  leading  edge.  In  the  combustion  section  a  large  separation  bubble  is 
established  behind  the  injector  that  thins  to  a  narrow  region  of  separated  flow  near  the 
wall  of  the  combustor.  In  contrast,  a  larger  nozzle  diameter  (1.6  inches)  produces  large 
separated  flow  regions  that  dominate  both  the  inlet  and  the  combustion  sections. 

Figure  10  shows  comparisons  between  internal  surface  pressure  measurements  and 
turbulent/detached  flow  computations  for  three  SFRJ  configurations  with  an  injector  di¬ 
ameter  of  1.9  in.  and  nozzle  diameters  of  1.6  in.,  1.54  in.,  and  1.1  in.  The  criterion 
established  in  the  previous  section  was  used  to  specify  the  flow  trainsition  location.  The 
similarity  of  these  results  with  those  for  the  1.7  in.  nozzle  (Figure  7)  indicates  the  relative 
insignificance  of  the  injector  diameter  in  determining  the  internal  wall  pressures. 


IV.  CONCLUSIONS 


Computational  solutions  of  the  internal  and  external  flowfield  for  inert  SFRJ  pro¬ 
jectiles  at  zero  yaw  have  been  performed  using  a  zonal  Navier-Stokes  code  with  a  TVD 
scheme.  Comparisons  between  computed  and  measured  internal  surface  pressures  were 
made. 

The  flow  transition  location  was  guided  by  measured  internal  wall  pressure  data  and 
chosen  to  achieve  the  best  overall  agreement  with  measured  data.  While  the  inlet  flow 
transition  point  is  determined  by  the  expelled  normal  shock  for  subsonic  internal  flow  (1.1 
inch  nozzle  diameter),  for  supersonic  flow  (1.6  inch  nozzle  diameter),  the  transition  point 
location  is  unclear.  The  best  agreement  with  measured  wall  pressures  is  achieved  when  flow 
transition  is  specified  at  the  end  of  the  inlet  (at  the  injector).  The  best  overall  agreement 
with  the  measured  pressure  level  and  gradient  along  the  combustion  section,  was  achieved 
using  the  backflow  turbulence  model  of  Goldberg. 

The  agreement  with  measured  internal  wall  pressures  lends  assurance  that  computed 
flow  predictions  can  correctly  replicate  physical  flow  details.  This  agreement  was  achieved 
with  the  aid  of  experimental  data  to  guide  the  location  of  flow  transition.  These  com¬ 
putations  indicate  that  nozzle  diameters  of  about  1.6  inches  yield  a  leading  edge  oblique 
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shock  and  supersonic  core  flow.  The  inlet  and  combustion  sections  are  dominated  by  flow 
detachment.  A  nozzle  diameter  of  1.1  inches  yields  an  expelled  leading  edge  normal  shock, 
subsonic  internal  flow,  and  a  choked  nozzle. 
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Figure  2.  Zone  Designations  for  Axisymmetric  SFRJ  Geometry.  Dimensions  in  Inches. 
1.  Inlet,  2.  Injector;  3:  Combustion  Section;  4:  Nozzle;  5:  Base;  6:  Exterior. 
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Na vier-Stokes  Computation  for  SFRJ 

1.9  in.  Injector,  1.1  in.  Nozzle  =  4.03 


Figure  6.  Internal  Surface  Pressure  Distribution  for  SFRJ,  il/oo=4.03,  1.9  in.  Injector, 
1.1  in.  Nozzle.  Effect  of  Flow  Transition  Location. 


Navier-Stokes  Computation  for  SFRJ 

1.7  in.  Injector  =  4.03 


7’  I"t"n&1  S"rface  Pressure  Distribution  for  SFRJ,  A/^4.03,  1.7  in.  Injector 
1.1,  1.54,  and  1.6  in.  Nozzles. 


-0.50  2.26  5.02  7.78  10.54 


Figure  8.  Pressure  Contours  for  SFRJ,  Moc=4.03.  1.7  in.  Injector,  1.6  in.  Nozzle. 
Dimensions  in  Inches. 
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Figure  9. 

Dimensions 


Pressure  Contours  for  SFRJ,  Mx= 4.03.  1.7  in.  Injector.  1.1  in. 
in  Inches. 


Nozzle. 
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Na vier-Stokes  Computation  for  SFRJ 

1.9  in.  Injector  M_  =  4.03 
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O  1.6  in.  Nozzle  W.T.  Data 
1.6  in.  Nozzle  Computed 
A  1.54  in.  Nozzle  W.T.  Data 
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Figure  10.  Internal  Surface  Pressure  Distribution  for  SFRJ,  Mx= 4.03,  1.9  in.  Injector, 
1.1,  1.54,  and  1.6  in.  Nozzles. 
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